Break-down of the single-active-electron approximation for one-photon ionization of 
the B X S+ state of H 2 exposed to intense laser fields 
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Ionization, excitation, and de-excitation to the ground state is studied theoretically for the first 
excited singlet state B X E+ of H2 exposed to intense laser fields with photon energies in between 
about 3eV and 13 eV. A parallel orientation of a linear polarized laser and the molecular axis is 
considered. Within the dipole and the fixed-nuclei approximations the time-dependent Schrodinger 
equation describing the electronic motion is solved in full dimensionality and compared to simpler 
models. A dramatic break-down of the single-active-electron approximation is found and explained 
to be due to the inadequate description of the final continuum states. 
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I. INTRODUCTION 

Recently, there has been an increased interest in the 
development of intense light sources with wavelengths 
shorter than the one of the titanium:sapphire lasers. This 
includes both high-harmonic sources and, especially, the 
free-electron lasers (FEL) [1]. Although the primary goal 
is to achieve radiation in the x-ray regime, also the vac- 
uum UV is of some practical interest. Interesting results 
have recently been achieved, e.g., at the free-electron 
laser FLASH in Hamburg [2-5]. 

For a laser pulse with a photon frequency centered 
at ijj and a peak electric field Fq the Keldysh parame- 
ter 7 = Ip/Up [6] with the ionization potential I p and 

F 2 

the pondcromotive energy U p — is usually used as a 
(rough) criterion to distinguish the so-called multiphoton 
(7 > 1) and tunneling (7 < 1) regimes. With the large 
photon energies as they are available from FELs one is 
likely to remain in the multiphoton regime. Interesting 
questions in this context are whether at these photon 
energies even lowest-order perturbation theory (LOPT) 
may be sufficient to properly describe most ionization 
processes in such laser pulses and at which intensities 
non-pcrturbative behavior sets in. A natural candidate 
for investigating these questions is H 2 , since its two elec- 
trons exposed to intense laser pulses can be treated per- 
turbatively [7, 8] and non-perturbatively [9-14] within 
full dimensionality. 

Despite the fact that electronically excited states have 
been found theoretically and experimentally to be popu- 
lated by a non-negligible amount in intense laser pulses 
[15-21], most theoretical studies on H 2 exposed to laser 
fields have concentrated on the ground state as initial 
state. However, the field-induced coupling to the typi- 
cally much more closely spaced neighbor states may com- 
plicate the strong-field behavior of excited states. On the 



other hand, since higher photon fluxes require lower tar- 
get densities, the novel intense light sources may allow 
for direct studies of electronically excited states similar as 
for ions. Such an investigation was, e. g., recently demon- 
strated for HeH+ [22, 23]. 

To the authors' knowledge, the only theoretical studies 
of excited states of H 2 in intense laser fields were reported 
in [24] in which their behavior in a (quasi)static field was 
investigated and in [25] in which one-photon ionization 
of the B 1 E+ state was considered within a newly devel- 
oped single-active-electron (SAE) approximation based 
on Koopmans' picture. In fact, the latter work has moti- 
vated the present investigation. Thus the first electroni- 
cally excited state B 1 !]^ of H 2 exposed to intense laser 
fields is considered and the validity of the SAE approx- 
imation for this excited state is investigated. The H 2 
molecule is treated in the fixed-nuclei approximation at 
R = 2.2998 ao with the field polarization being parallel 
to the internuclear axis. First, 8eV photons were consid- 
ered as in [25]. Motivated by the found order of magni- 
tude deviation between the present full two-electron cal- 
culation and the one in [25] , a larger photon-energy range 
within the regime of one-photon ionization is considered, 
but only for the case of a single open channel, i. e., before 
the photon energy is sufficient to leave the H^ ion in an 
electronically excited state. It may be noted that some of 
the photon energies used in the calculations are already 
available with FELs [26, 27]. 

After a brief description of the methods in Sec. II the 
results are presented and discussed in Sec. III. This in- 
cludes a discussion of the results for a photon energy 
of 8 eV in Sec. Ill A and for variable photon energies in 
Sec. Ill B. A simple model for explaining the failure of 
the SAE is given in Sec. Ill C, followed by a conclusion in 
Sec. IV. Atomic units are used, if not stated otherwise. 



II. METHOD 
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Most results of this work have been obtained by a full- 
dimensional solution of the time-dependent Schrodinger 
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equation (TDSE) describing the two electrons of H2 ex- 
posed to an intense laser pulse with linear polarization. 
The molecular axis is assumed to be: aligned parallel to 
the electric field of the laser. The semi-classical non- 
relativistic dipole approximation is adopted in which the 
laser field is described classically, while the molecular 
system is treated quantum mechanically. The shown 
results are obtained for fixed internuclear separations. 
The details of the approach have been given previously 
[10, 11, 13, 20] and are thus only briefly repeated. 

The electronic TDSE describing H 2 in a laser field is 
given within the above-mentioned approximations as 

»^*(n,r 2 ,t) = (Ho + D(t)) *(ri,r 2 ,t) (1) 

where Ho is the field-free electronic Hamiltonian of H 2 
and D(i) describes the interaction with the field. This 
interaction may be given in either velocity form, T)(t) = 
A(£)-(pi+p 2 ), or in length form, D(t) = — F(i)-(r 1 + r 2 ), 
where Pi is the momentum operator of electron i while 
A and F represent vector potential and the electric field, 
respectively. 

The wavefunction 

*(ri,r 2 ,i) = J2 KL{t)4>nL{vi,v 2 ) (2) 

nL 

is expanded in terms of the field-free states <p n L and the 
time-dependent coefficients b n L(t). The field-free wave- 
functions 4>nL (ri,r 2 ) and corresponding energy eigenval- 
ues E n L are obtained from the solution of the time- 
independent Schrodinger equation (TISE), 

Ho0„i(ri,r 2 ) = E nL <t> nL {T U T 2 ) • (3) 

A combined index L is used for the specification of the 
symmetry. In this work it is limited to either 1 Ti g (L = 0) 
or 1 E tl (L — 1) symmetry. The index n — 1,2,... just 
numbers the different states for a given symmetry L. 

The field-free two-electron wavefunctions <p n L are ob- 
tained by a configuration-interaction (CI) calculation 
performed in a basis of orbitals. The orbitals 
are expressed in a B-spline basis set in prolate-spheroidal 
coordinates. Since the basis functions are confined to a 
finite spatial volume, a discretized representation of the 
electronic continuum is obtained. Substitution of Eq. (2) 
into Eq. (1) yields then a finite set of coupled differential 
equations for the time-dependent coefficients These 
equations are solved numerically with the initial condi- 
tion b n L(t = 0) = <5„,i Sl,i for the B 1 S+ state. 

In the present work, the TISE for is solved in a 
box of size 350 do- The solution of the TISE for is 
obtained in prolate spheroidal coordinates (1 < £ < 
00,— 1 < rj < 1,0 < <p < 2tt). The total symmetry of 
the states is defined on the basis of angular momen- 
tum and gerade or ungerade symmetry of the state wave- 
function. For each total symmetry of 350 B splines 
of order 15 are used in £ direction and 24 B splines of 



order 8 in i] direction. H^" states with angular momenta 
between and 5 with gerade and ungerade symmetry are 
calculated. Using such a basis set yields around 4200 
states for each total symmetry of H \ . The CI calcula- 
tion performed with these states gives around 8000 
states for each symmetry of H 2 and 

For computational convenience, cos 2 -shaped pulses 

X(t) = X cos 2 (nt/r) cos(wi + C ep) (4) 

where lj is the photon energy and X stands for F (length 
form) or A (velocity form) are used in most of the calcu- 
lations of this work. (Note, the pulses defined this way 
are not identical for X — F or X — A, especially in the 
case of extremely short pulses.) The peak electric field or 
vector potential amplitudes are F and A , respectively. 
In order to obtain the maximum field at the center of the 
pulse, the carrier-envelope phase (CEP) </>cep is set to 
(71-/2) for the length (velocity) form. The advantage of 
the cos 2 pulses is the finite pulse length r (in the interval 
[— r/2, t/2]) and thus the well-defined interval for time 
integration. This is an evident advantage compared with 
the (more realistic) Gaussian pulse 

( 2 i 2 \ 

X(t) = X cxp y-—J cos(uit + <Pcep) (5) 

that only exponentially decays to zero and requires thus 
a careful convergence study with respect to the integra- 
tion interval. In order to assure that found deviations to 
the calculation in [25] are not caused by a possible dif- 
ference in the definition of the pulse shape, the results 
of cos 2 and Gaussian pulses are compared. In this case 
the characteristic time r g of the Gaussian pulse is chosen 
according to 

This choice yields the same full width at half maximum 
(FWHM) of the Gaussian and the cos 2 pulses. 

For sufficiently low intensities one expects lowest-order 
perturbation theory (LOPT) to adequately describe the 
interaction of a molecule with a laser field. Within LOPT 
the TV-photon ionization rate (in s _1 ) is given by 

r<») . „ m (£f (7) 

where / is the intensity of the monochromatic laser field 
in W/cm 2 , hw is the photon energy in Joule, and 
is the generalized TV-photon ionization cross section in 
cm 2Ar s A ' _1 . For single-photon ionization is equal to 
the standard one-photon cross-section and LOPT reduces 
to Fermi's Golden Rule. Integration of Eq. (7) over the 
intensity profile of a pulse gives the ionization yield 

P^ n OPT = 1 - exp (- f T^(I(t))dt) (8) 

V JPulsc / 
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where is a function of intensity I (see Eq. (7)) which 
for a laser pulse is a function of time t. The generalized 
cross-sections are in this work calculated with the 
same field-free two-electron states 0„l of Eq. (3) that are 
used for solving the TDSE. Results within LOPT have 
been earlier obtained for H2 in [7, 8] and a comparison 
to TDSE results was given in [10]. In those references 
technical details on the evaluation of a^ N>> may be found. 
However, in those previous investigations only ionization 
from the electronic ground state was considered. 

For a more detailed investigation of the validity of the 
SAE approach introduced in [25] it is of interest to con- 
sider also alternative SAE models. In one approach the 
TDSE is solved as discussed above, but a restricted set of 
configurations is used in the CI calculation. In this case 
only configurations are adopted in which one electron oc- 
cupies the lowest lying (la g ) orbital of H^. The other 
electron may then occupy any of the na g (na u ) orbitals 
to yield the configurations for the CI calculation giving 
the ( £„) states of H2. This approach is referred 
to as pseudo SAE (p-SAE) [19, 20, 28]. The p-SAE ap- 
proach corresponds to a complete relaxation of the core. 
A frozen-core SAE approach where there is no relax- 
ation of the core has recently been introduced in [20]. 
In that case the core electrons are frozen in their field- 
free initial-state orbitals and the TDSE is solved for the 
active electron in the combined field of the core and the 
laser pulse. In the implementation the core was described 
within density- functional theory (DFT). For strong- field 
excitation and ionization of ground-state H 2 it was, how- 
ever, demonstrated that the results agree quite well with 
the ones obtained with a core described within Hartree- 
Fock theory, despite the rather different electronic bind- 
ing energies. The DFT-based variant used in this work 
will be referred to as DFT-SAE. 

Two of the SAE approaches, p-SAE and DFT-SAE, 
give the ionization yield for a single electron. In [20] it 
was found that for ionization yields of up to about 10 % 
the SAE results for H 2 should be multiplied by a factor 2 
in order to properly account for the two equivalent elec- 
trons. However, for larger ion yields this factor 2 leads 
to an overestimation, because the screening of the core 
electron is reduced. As a consequence, the ionization po- 
tential increases and the ionization probability decreases. 
Therefore, the ion yield of the SAE approximation ap- 
proaches the full CI-TDSE result for high intensities and 
ion yields larger than about 20%. However, this screen- 
ing argument applies in principle only, if the ionization 
process depends strongly on the electronic binding en- 
ergy. This is, e.g., the case for ionization in the tunnel- 
ing picture, but is not expected to be valid in the op- 
posite extreme of perturbative single-photon ionization. 
It should be noted that no prefactor should be used in 
the Koopmans'-picture based SAE approach (K-SAE) of 
Barth et al. [25]. 



III. RESULTS 

In Table I the energies of the ground (X 1 ^!") and the 
first excited (B £+) states of H 2 are given at the in- 
ternuclear separation R — 2.2998a . Furthermore, the 
corresponding ground-state energies of the Hj ion are 
shown together with the ionization potential of the B 
state that follows from them. The energies obtained by 
means of K-SAE (given in [25]), CI, p-SAE, and DFT- 
SAE approaches are compared. The best available the- 
oretical energy values are also shown. While the energy 
of the B *£+ state obtained with K-SAE is m very good 
agreement with the CI calculation, the DFT-SAE result 
is almost 0.09 a.u. (2.4 eV) away. The p-SAE result is 
much closer to the CI result, but not as close as K-SAE. 
Since the ionic ground-state energies agree very well for 
all approaches, the differences found for the ionization 
potentials are a result of the different B : S+ energies. 
For completeness, also the different ground-state ener- 
gies are given, if available. In this case, DFT-SAE agrees 
better to CI than p-SAE. It was demonstrated in [31, 32] 
that very accurate CI results can be obtained with the 
present approach, if the basis-set parameters including 
the configurations are judiciously chosen. However, for 
solving the TDSE a compromise has to be made, since a 
large spectrum of field-free states including the electronic 
continuum has to be described with reasonable accuracy. 

The reason for the very poor DFT energy of the B 1 £+ 
state is probably twofold. First, small systems like He 
and H 2 are known to be difficult to be described by DFT. 
More importantly, however, standard DFT is in princi- 
ple not applicable to excited states. In practice, time- 
dependent DFT is usually adopted to obtain excitation 
energies. Nevertheless, it should be interesting to inves- 
tigate the performance of the recently developed DFT 
based SAE in comparison to the also very newly pro- 
posed K-SAE as is done below. 

A. 8 eV photons 

In [25] it was concluded that single-photon ionization 
of the B 1 £+ state of H 2 at 8 eV should satisfy the condi- 
tions for the applicability of a time-dependent extension 
of Koopmans' picture, a single-active-electron approxi- 
mation based on Koopmans' theorem, i. e. within K-SAE. 
Therefore, the response to a 5-cycle laser pulse with 8eV 
photons was studied in the framework of the newly devel- 
oped K-SAE. The TDSE was solved on a grid using the 
length formulation. In [25] the calculated ionization yield 
plotted against the laser peak intensity on a log-log scale 
showed a slope 1 in a large range of intensities. This con- 
firmed the perturbative nature of the single-photon pro- 
cess in the selected intensity window and was interpreted 
as a further confirmation of the proper implementation of 
the time-propagation in the K-SAE approach. Note, in 
[25] the laser peak intensity is defined as / = cqcFq where 
eq is the vacuum permittivity, c is the vacuum speed of 
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Table I. State energies (in atomic units E h ) using different approaches at R — 2.2998 a . The ionization potential of the B 1 E+ 
states is shown in the last column. 

Approach H 2 X 2 S+ H 2 H 2 + X% E^B 1 ^) 

CI -1.10374 -0.7485 -0.5989 0.1496 

p-SAE -1.06385 -0.7388 -0.5989 0.1399 

DFT-SAE -1.08196 -0.8381 -0.5989 0.2392 

K-SAE [25] n.a. -0.7509 -0.5988 0.1521 

Accurate -1.1111725 [29] -0.7563608186 [30] -0.5989 0.1575 



light, and F is the peak value of the electric field. In 
this work the alternative, cycle averaged peak intensity 
/ = ^e c\F 2 is used. 

Table II. Comparison of K-SAE, CI, LOPT, p-SAE, and DFT- 
SAE results. The populations of the initial B 1 Ej state and 
the ionization yields are given for two different values of the 
peak electric field. The K-SAE results are taken from [25]. 





F (a.u.) 


K-SAE 


CI 


LOPT 


p-SAE 


DFT-SAE 


Initial 


0.005 


n.a. 


0.960 


0.963 


0.967 


0.980 


Ion. 


0.005 


0.0043 


0.033 


0.037 


0.0031 


0.022 


Initial 


0.020 


0.930 


0.548 


0.543 


0.624 


0.728 


Ion. 


0.020 


0.066 


0.379 


0.457 


0.036 


0.261 



Besides the graphically shown intensity scan, the ion- 
ization yield and (in one case) the population remain- 
ing in the initial B 1 !]^ state of H2 were given in nu- 
merical form for the peak electric fields F = 0.005 and 
0.02 a. u. in [25]. In Table II these results are compared 
to the corresponding values obtained in this work within 
the different approaches, CI-TDSE, LOPT, p-SAE, and 
DFT-SAE. In LOPT the initial-state population P in i was 
defined as P ini = 1 - i^° PT 

The ionization yield predicted on the basis of K-SAE 
differs substantially from the CI-TDSE result. At the 
lower field strength (F = 0.005 a.u.) the K-SAE yield 
is almost one order of magnitude (factor 7.7) too small. 
For Fq — 0.02 a. u. the disagreement is a little bit smaller 
(factor 5.7). While p-SAE deviates from the CI-TDSE 
results even more than K-SAE (underestimation by fac- 
tors 10.6 and 10.5), DFT-SAE appears on the first glance 
to be the SAE model with the best agreement to the 
full two-electron result, despite the poor energy of the 
B 1 S+ state. Although the DFT-SAE yields are also 
smaller than the CI-TDSE ones, they differ only by fac- 
tors 1.5 and 1.45 for the smaller and the larger field 
strengths, respectively. Clearly the best agreement to 
the CI-TDSE result is, however, obtained with LOPT. 
At the lower field strength the agreement is in fact good 
(the CI-TDSE yield is overestimated by a factor 1.1 
only) and for the larger field strength it is still reason- 
able (factor 1.2). Note, the given fields correspond to 



laser peak intensities of about 8.8 x 10 11 W/cm 2 and 
1.4 x 10 13 W/cm 2 where LOPT is not necessarily ex- 
pected to work. The one-photon cross-section at 8cV 
is found to be crW = 5.75 x 10" 17 cm 2 . 

A Taylor expansion of the exponential in Eq. (8) and 
insertion of Eq. (7) shows that for one-photon ionization 
the ion yield is proportional to the one-photon absorption 
cross-section and the laser intensity, f^O PT (x pi 1 ) j if 
the time integral is small. Therefore, one expects within 
LOPT and for small ion yields that an increase of the 
field strength by a factor 4 leads to an increase of the ion 
yield by a factor 4 2 = 16. This factor is reproduced well 
for K-SAE (15.3), but neither for CI-TDSE (11.5) nor for 
LOPT (12.4). The latter result shows that for the larger 
ionization yields found for the two-electron models CI- 
TDSE and LOPT the assumption of a sufficiently small 
time integral is not fulfilled and the ion yield does not 
increase linearly with intensity, even within LOPT. Cor- 
recting the K-SAE yield at the larger field for this effect 
shows that the agreement between K-SAE and CI-TDSE 
does in fact not improve for the higher field, but the dis- 
agreement remains at a factor of about 7.5. While the 
factor 11.9 found between the two ion yields obtained 
with DFT-SAE appears understandable from the simi- 
lar magnitude of ionization as compared with CI-TDSE 
and LOPT, one would expect for p-SAE a factor close to 
16 similar to K-SAE, since the ion yield is even smaller. 
However, instead a value 11.6 is found, similar to the 
other cases except K-SAE. 

The pronounced failure of K-SAE despite the very ac- 
curate energy of the B 1 E+ state is clearly surprising. In 
order to clarify this issue it was checked that it is not an 
artifact of the different numerical implementations. One 
possibility could, e.g., be a different pulse definition, or 
the use of the velocity or length forms of the dipole op- 
erator. In Fig. 1 the ionization yield as a function of the 
laser peak intensity is shown. It confirms the already 
discussed good agreement of CI-TDSE with LOPT. The 
ionization yield obtained by means of DFT-SAE is quite 
close to the full CI-TDSE results. The p-SAE and K-SAE 
fail by almost an order of magnitude for 8eV photons. 
In order to investigate the sensitivity of the results on 
pulse details and the form of the dipole operator used, 
two CI-TDSE results are shown in Fig. 1. In one case 
a cos 2 pulse envelope (for E(t)) is used with dipole mo- 
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Figure 1. Ionization yield as a function of the peak laser inten- 
sity. The TDSE calculation is performed using length-form 
transition dipole moments with a cos 2 envelope for the electric 
field (black, circles) and velocity-form transition dipole mo- 
ments with a Gaussian envelope for the vector potential (red, 
squares). The FWHM for both pulses is 2.5 cycles of 8eV 
photon energy. The ionization yields obtained with LOPT 
(blue, solid line), DFT-SAE (violet, diamonds), scaled DFT- 
SAE by shifting the photon energy by — 2.44eV (at 10.44 eV) 
(orange, asterix), and p-SAE (maroon, plus) are also shown 
and compared to the results from Barth et al. [25] (green, 
triangles). 

merits in length form, while in the other case a Gaussian 
pulse shape (for A(t)) is used together with the velocity 
form. Clearly, neither the choice of the dipole operator 
nor the exact pulse shape (in both cases a FWHM of 2.5 
cycles is used) modifies the result substantially. In fact, 
the agreement is almost perfect, except for very high in- 
tensities and thus very close to saturation. 

In view of the rather poor energy of the B 1 S+ state in 
DFT-SAE the superior result compared to the other SAE 
models is suspicious. In fact, as has been discussed pre- 
viously in [8, 10, 33] a poor initial-state energy should 
be accounted for by a correspondingly shifted photon 
energy. This energy shift should be chosen in such a 
way that the correct multi-photon ionization threshold 
(in this case the one-photon threshold) is obtained. Since 
the present DFT threshold lies about 2.44 eV higher than 
the one of the CI calculation (see Table I), the CI-TDSE 



results for a photon energy of 8 eV should be compared 
to the DFT-SAE results at a photon energy of 10.44 eV. 
The corresponding curve (also drawn in Fig. 1) shows 
that the agreement with CI-TDSE worsens evidently, al- 
though it remains better than the agreement of the other 
SAE models with the CI-TDSE result. It is clear, how- 
ever, that the (slightly) better agreement of DFT-SAE 
appears to be accidental. 

It is, of course, interesting to investigate the origin 
of the rather evident failures of K-SAE and p-SAE. A 
characteristic property of the B 1 S+ state of H 2 is its 
rather large ionic component [34]. In fact, it was shown 
theoretically that this ionic character plays a crucial role 
in the strong-field behaviour of HV For example, it is 
responsible for bond softening and enhanced ionization 
of neutral H2 [9, 35-38]. The field- induced admixture of 
the ionic (H + H~) component to the X 1 E+ ground state 
leads to strongly enhanced ionization, since H _ possesses 
a very low electronic binding energy. It appears quite 
reasonable that SAE models fail to properly describe the 
ionic component of the B 1 E+ state, since the (frozen) 
spectator electron always occupies a molecular orbital 
that is symmetrically distributed over both protons. As 
already discussed, the lack of even a very small admixture 
of the ionic component may reduce the ionization rate in 
a strong electric field significantly. This could explain 
the too low ionization yield found uniformly for all SAE 
models discussed in this work. 

However, the good agreement between LOPT and CI- 
TDSE (especially at low intensities) allows to conclude 
that one should not look for a strong-field explanation 
where the binding energy influences the ionization rate 
in an exponential way. Instead, the validity of the LOPT 
indicates that it is only the (direct) transition dipole ma- 
trix element between the initial B 1 E+ state and the con- 
tinuum (reached with an 8 eV photon) that determines 
the magnitude of ionization. In fact, a rough estimate 
based on a simplified model describing ionization from 
such an ionic pair state (weighted by the contribution of 
this state to the B X E+ state of H2) indicates that this 
ionic component cannot at all explain a massive increase 
of the one-photon ionization yield. Another possible ex- 
planation of the failure of a SAE model could be the 
occurrence of doubly-excited autoionizing states that are 
necessarily absent in SAE models. However, they are 
expected to be excited at photon energies larger than 
8eV[32, 39]. 

B. One-photon ionization 

A more complete picture of one-photon ionization from 
the B 1 E+ state of H2 may be obtained from a consider- 
ation of the photon-energy dependence. Figure 2 shows 
the ionization yield of the H2 B 1 E+ state exposed to a 
10 fs laser pulse with a peak intensity of 1.0 x 10 9 W/cm 2 
obtained with CI-TDSE. The photon energy varies be- 
tween about 3 and 13 eV and covers thus the one-photon 
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Figure 2. Ionization (black, circles), population of the other 
bound states OBS (red, dashed), and population of the X x Eg 
ground state (blue, asterix) obtained with CI-TDSE as a func- 
tion of the incident photon energy for a 10 fs laser pulse with 
a peak intensity of 1.0 x 10 9 W/cm 2 . 



regime from just below the ionization threshold to the 
energy at which the first excited electronic state of 
becomes energetically accessible. Due to the finite pulse 
width (and possible two-photon ionization) the ioniza- 
tion starts before the threshold energy of about 4.1 eV. In 
contrast to one-photon ionization starting from the elec- 
tronic ground state which shows an overall monotonous 
decrease with increasing photon energy [40], the ioniza- 
tion yield of the B 1 E+ state decreases until about 5eV 
after which it increases drastically (note the logarithmic 
scale) until about 7 eV. Beyond about 7 eV the yield de- 
creases monotonously, until some structure due to the 
occurrence of autoionizing states becomes visible in be- 
tween about 11 eV and 12.5 eV. Clearly, the autoioinizing 
states cannot explain the failure of the SAE approaches 
at 8 eV discussed in Sec. Ill A. 

The population of the H2 X 1 S+ ground state after the 
laser pulse is also shown in Figure 2. A pronounced max- 
imum occurs at the expected energy of about 9.66 eV at 
which the B 1 E+ — > X 1 E+ transition becomes resonant. 
The width of the peak reflects the spectral width of the 
laser pulse. For the considered laser peak intensity the 
population of the ground state and thus the probability 
for deexcitation of the B 1 ^^ state is at its maximum 
about a factor 10 larger than the ionization yield. Since 
the total depopulation of the B 1 !]^ state is, however, 
still very small (less than 0.1 %), the ionization probabil- 
ity remains practically uninfluenced by this deexcitation 
process. As was shown recently in [13], the oscillations 
(side peaks) of the resonant peak are due to the cos 2 
shape of the pulse used in the calculation. Since they 
are many orders of magnitude smaller than the ioniza- 
tion yield, the latter is, however, not influenced by these 
spurious oscillations. This is also confirmed explicitly for 
8eV photons in Fig. 1, since the ion yields obtained with 
a Gaussian pulse that does not show these oscillations 
agrees well with the ones obtained with a cos 2 pulse. 



The population of the other bound states (OBS) is also 
shown in Figure 2. It is denned as Pobs = 1 — (Pni+Pon) 
where p n i is the population of the initial B state and 
Pi on is the ionization yield). In a large range of photon 
energies Pobs is dominated by the ground-state contri- 
bution. Only below 6 eV the energetically higher lying 
bound states come into play. Close to the ionization 
threshold the ion yield and Pobs cross each other. Below 
threshold the Rydberg states of H2 are resonantly popu- 
lated, while the one-photon ionization channel closes. As 
already mentioned, the finite band width of the adopted 
laser pulse leads to a rather broad energy range in which 
excitation and ionization compete with each other. 




6 8 10 

Photon energy (eV) 

Figure 3. Photon-energy dependent ionization yield for var- 
ious theoretical models with laser parameters as in Fig. 2. 
The ionization yields obtained with CI-TDSE (black, circles), 
LOPT(red, dashed), p-SAE (blue, square), and DFT-SAE 
(violet, diamonds) are compared. The ionization yields ob- 
tained by a shift of the photon energy to compensate the 
corresponding error in the ionization potential with respect 
to CI of p-SAE (maroon, triangles) and DFT-SAE (green, 
plus) are also shown. 

Figure 3 shows the ionization yields obtained by means 
of the different theoretical models in the one-photon ion- 
ization regime. The laser parameters are the same as 
the ones in Figure 2. The ionization yield obtained by 
CI-TDSE as shown in Figure 2 is repeated for compari- 
son. For the chosen laser intensity the Keldysh parameter 
is much larger than 1 even at the lowest photon energy. 
This indicates ionization to take place in the multiphoton 
regime (although the name is evidently a bit misleading 
in the present case of one-photon ionization). In fact, 
as the comparison to the prediction of LOPT shows, the 
process is quite well described by perturbation theory in 
almost the complete energy range. The observed devia- 
tions close to the ionization threshold and in the regime 
of autoionizing states are mainly due to the finite pulse 
width in the TDSE calculation. Since LOPT predicts 
ionization rates, an infinite pulse duration is implied. As 
a consequence, a sharp ionization threshold and sharp 
resonances with the widths of the latter determined only 
by the corresponding lifetimes are observed in the LOPT 
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spectrum. A more appropriate comparison would thus 
involve to convolute the LOPT spectrum with the spec- 
tral band width or to consider longer pulses in the TDSE 
calculation. In fact, it is interesting to note that already 
for a pulse as short as 10 fs such a good agreement is 
found. For the present discussion the most important 
conclusion is, however, the very good agreement of LOPT 
and CI-TDSE in the non-resonant photon-energy range 
between about 6 and 10.5 eV in which sharp features are 
absent. 

The ionization yield obtained by solving the TDSE 
in the p-SAE approximation shows beyond the ioniza- 
tion threshold a monotonous decreasing behavior with 
increasing photon energy. This is very similar to the 
behavior found for one-photon ionization from the elec- 
tronic ground state. In order to correct for the inaccurate 
ionization potential obtained within the p-SAE approach 
the ionization yield may be shifted by 0.264 eV in order 
to compensate this effect. However, even after the shift 
the ion yields obtained with the p-SAE-TDSE approach 
or CI-TDSE differ substantially. At the threshold the p- 
SAE result is about one order of magnitude larger than 
the CI yield. At about 6 eV the ion yields of the p-SAE 
and CI approaches cross each other. At the maximum 
of the CI ionization yield close to 7.3 eV the CI result 
is about an order of magnitude large than the p-SAE 
result. Apart from the local variation in the resonant 
regime, the ratio between CI and p-SAE results decreases 
then rather monotonously to about a factor 2.5 at 13 eV. 
Clearly, p-SAE is completely inadequate for predicting 
the one-photon ionization yield of the B state of H2 
in the whole considered photon-energy range, and not 
only at 8 eV. 

The DFT-SAE ionization yield shifted by -2.44 eV 
in order to compensate the wrong ionization potential 
agrees at the threshold, i.e. between about 3.5 eV and 
4.75 eV quite well with the p-SAE result. However, al- 
though the DFT curve decreases also monotonously with 
increasing photon energy, its decrease is much slower 
compared to the one predicted within the p-SAE ap- 
proach. As a consequence of the smaller slope, the quan- 
titative deviation between DFT-SAE and the CI calcu- 
lation is smaller than the one between p-SAE and CI. 
In view of the massive qualitative difference between the 
DFT-SAE and the CI curve this could, however, be ac- 
cidental. Clearly, the reasonable agreement between the 
unshifted DFT-SAE and the CI result at 8eV is defi- 
nitely accidental. In conclusion, Fig. 3 shows that both 
SAE approaches (p-SAE and DFT-SAE) fail completely 
when describing one-photon ionization of the B 1 E+ state 
of H2 . Quantitatively, there is a disagreement of up to 
an order of magnitude. Due to the completely differ- 
ent qualitative behavior, the SAE results may under- or 
overestimate the true ionization yield. The result of the 
K-SAE at 8 eV [25] seems to indicate, that a similar fail- 
ure as is found for p-SAE is to be expected also for the 
K-SAE, despite the relatively accurate initial-state de- 
scription of the latter. In fact, as is shown in the next 



section, it is really a failure of the SAE approaches to 
properly describe the final continuum states. 

C. Explanation for the failure of the SAE 
approximation 




Photon energy (eV) 



Figure 4. The one-photon ionization cross-section (within 
LOPT) of the H 2 B 1 ^ state for the CI (solid, black), the 
CP (red, dashed), and the p-SAE (blue, chain) basis sets. 

In order to analyze the origin of the massive devia- 
tion between the SAE and the CI ionization yields it 
may be noted that the p-SAE approach differs from the 
full CI calculation only by a restriction of the configu- 
rations that are included in the CI calculation, as was 
outlined in Sec. II. Furthermore, the previous results 
have demonstrated that LOPT is applicable and the de- 
viations between the CI and the p-SAE LOPT cross sec- 
tions is practically identical to the ones found between 
CI-TDSE and p-SAE-TDSE. However, LOPT allows a 
much easier analysis of the failure of the SAE approxi- 
mation, as is outlined below. 

The origin of the failure of the SAE approaches can be 
understood from Fig. 4 that shows the one-photon ioniza- 
tion cross-sections within LOPT obtained from CI calcu- 
lations with 3 different basis sets. The CI basis that gave 
converged results is again denoted as CI. As already ex- 
plained, keeping in the CI basis only those configurations 
in which one electron remains in the 1 a g ground-state or- 
bital of Hj leads to the p-SAE model. The CP basis set 
is almost identical to the p-SAE basis set. The only dif- 
ference is that it includes one additional |1<7 U )|1<7 U ) con- 
figuration in the CI calculation for the states with 1 E+ 
symmetry. Thus the initial B 1 !]^ state is identical for 
the CP and the p-SAE basis sets. Only the final con- 
tinuum states differ by the additional admixture of the 
single configuration (|lcr u )|lcr u )). 

As is seen from Fig. 4, the addition of this configura- 
tion to the CI calculation of the final states leads to a 
substantial modification of the ionization cross-section. 
While for photon energies below about 6.5 eV the cross- 
section obtained with the CP is reduced compared to 
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the p-SAE result, there is a pronounced enhancement by 
more than an order of magnitude at higher photon ener- 
gies. While there is some evident deviation between the 
full CI calculation and the CP results, it is clear that the 
main failure of the SAE approximation is due to the lack 
of the |1<7„)|1<7 U ) configuration in the description of the 
final continuum states. A further inclusion of configura- 
tions is required for achieving quantitative convergence, 
but the cross section does not change qualitatively. The 
reduction at small energies becomes more pronounced 
and both minima and maxima are shifted to lower pho- 
ton energies. Clearly, also the adequate description of the 
spectral features above lleV due to autoionizing states 
requires the inclusion of additional configurations. 

Based on this finding a simple semi-quantitative model 
for the failure of the SAE approximation can be devel- 
oped. The one-photon ionization cross-section is pro- 
portional to |(?/>/ 1 D'li/'i)! 2 , where tpf and ipi are the final 

and the initial state wavcfunctions, respectively. D is 
in this case the time-independent electronic dipole mo- 
ment operator. In the simplest SAE approximation in 
which one electron is frozen during the whole process 
the initial state B^^S^ can be described as a product 
of the (mean-field) wavefunctions \l<r g ) and \la u ) of H2. 
A more sophisticated SAE model accounts for the proper 
(anti-) symmetry of the electronic wavefunction (as in the 
p-SAE model), but this does not change any of the fol- 
lowing conclusions as can easily be verified. In order to 
keep the notation simple, the simple frozen-core model 
is thus pursued. The final state ipf can then be approx- 
imated as \l(J g ) \eu g ) in this simple SAE approximation 
and as Ci(|la 9 ) |e<7 9 )) + C2{\l<J u )\la u )) in a simplified 
two-electron model with (some) correlation. Here, ea g 
describes the corresponding continuum orbital with en- 
ergy e and one has |Ci| 2 + |C 2 1 2 = 1. The dipole tran- 
sition matrix element in the SAE case is 

(V/|D|Vi> = (1<7 S ; ^ 9 |D|1<t 9 ; la u ) = (ea g \B\la u ) (9) 

while the two-electron model yields 

(V/IDIVO = Ci(«7 g |D|l<7„) + C 2 (la„|D|la s ) . (10) 

Evidently, the latter gets an additional contribution from 
the doubly excited \la u )\l<r u ) configuration. 

Of course, the size of the difference between Eqs. (9) 
and (10) depends on the magnitude of C2 and the rel- 
ative difference between (e<j g \t)\\(j u ) and {\a u \t)\la g ) . 
Although correlation in the final state and thus C2 is 
expected to be small, it can be compensated by the size; 
of the dipole matrix element. Using the additional sim- 
plification of adopting H J orbitals instead of the ones of 
H2 one finds (with the CP basis) the maximum value 
of the coefficient C2 to be about 0.15, and thus the 
maximum contribution of the | 1<t u ) | 1<t u ) configuration to 
the continuum wavefunction to amount to about 2.25%. 
While this coefficient is very small, the dipole moment 
(1(7^ |D| l<7g) between two bound orbitals of similar spa- 
tial extension is on the other hand much larger than the 



one between the continuum orbital ea g and the bound 
state la u . While the former is found to be equal to 
1.1602 a. u. at R = 2.3 ao (in good agreement with the 
results in [41, 42]), the discretized (and thus equally nor- 
malized) transition dipole matrix element (na g \D\la u ) 
is found to be 0.0408 a. u. at a final-state energy corre- 
sponding to a photon energy of about 8 eV. This explains 
easily the dominance of the second amplitude in Eq. (10) 
and the failure of an SAE model. 

Even a very small admixture of states like |lcr u )|l(7„) 
due to correlation can substantially modify the one- 
photon cross section, since the bound-bound transi- 
tion matrix elements are much larger than the bound- 
continuum ones. From Fig. 4 one may further conclude 
that the coefficient C2 changes sign (relative to C\) at 
a final-state energy at about 2 to 3 eV above the ioniza- 
tion threshold, since the cross section is first reduced and 
than enhanced compared to the one obtained with the 
SAE model. Clearly, a comparable effect does not occur 
for the electronic ground state, since all matrix elements 
of the type (n<r g ; mcr u \i>\la g ; la g ) vanish systematically 
for n ^ 1 due to the orthogonality of the orbitals; the 
matrix element of a one-electron operator between two 
Slater determinants differing by more than one orbital is 
zero. 



IV. CONCLUSIONS 

On the basis of results obtained from the solution 
of full dimensional time-dependent Schrodinger equation 
describing the two electrons of H 2 exposed to an intense 
laser field the ionization of the first electronically excited 
B 1 S+ state is investigated in the one-photon regime. A 
large deviation is found compared to a corresponding re- 
cent calculation [25] performed within a Koopmans' pic- 
ture based SAE approximation. It is demonstrated that 
also two other SAE approximations (called p-SAE and 
DFT-SAE) fail similarly, although DFT-SAE seemed on 
the first glance to yield reasonable results. However, this 
turned out to be accidentally caused by two partially 
compensating errors. On the other hand, lowest-order 
perturbation theory was shown to be quite adequate to 
reproduce the full TDSE calculations. This allowed to 
demonstrate that the failure of the SAE approximations 
is solely due to an improper description of the final con- 
tinuum state, while the initial B 1 !^ state is properly 
described. The inclusion of a single doubly-excited con- 
figuration into the configuration-interaction calculation 
of the final states explains semi-quantitatively the in- 
crease of the cross section by about an order of mag- 
nitude compared to an SAE model. The present work is 
thus an interesting example for a dramatic effect of cor- 
relation in the final continuum states that emerges for 
an excited initial electronic state but is absent, if the 
initial state is the electronic ground state of H 2 . The 
found effect may be observable with synchrotron radia- 
tion. However, since it is rather difficult to achieve the 
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required target density of electronically excited states, 
experiments with higher photon fluxes as they are avail- 
able from free-electron lasers should be advantageous. 
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